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SECTION  I 


INTRODUCTION 

While  there  has  been  great  progress  In  recent  years  in  the 
calculation  of  unsteady  transonic  flows rarely  have  any  of  the 
computational  procedures  which  have  evolved  been  applied  to  a set 
of  similar  problems  for  the  purpose  of  analytical  comparison.  The 
objective  of  this  study  is  to  provide  such  a comparison  for  two  of 
these  methods  in  order  that  the  character  of  the  resultant  solutions 
for  each  may  be  examined.  As  the  ultimate  use  of  any  flow  field 
computational  procedure  is  to  establish  the  aerodynamic  forces  nec- 
essary for  performance  evaluation  and  structural  analysis,  it  is  the 
intention  here  to  indicate  the  expected  behavior  of  calculations  for 
such  applications. 

Although  many  unsteady  transonic  flow  problems  have  been  con- 
sidered to  date,  there  exist  only  three  approaches  by  which  unsteady 
aerodynamic  forces  may  be  obtained;  namely,  harmonic  analysis,  time 
integration,  and  the  indiclal  method.  The  Indiclal  method,  while 
having  certain  inherent  advantages  for  structural  applications,  re- 
quires the  flow  field  response  to  a step  change  in  a given  mode  of 

(2) 

motion  of  the  airfoil  or  body  under  consideration  . This  "in- 
diclal response"  is  obtained  via  time  integration,  and  thus  this 
method,  in  spirit  at  least,  will  not  be  considered  distinct.  We, 
therefore,  limit  our  attention  to  a comparison  of  the  remaining  two 
procedures.  Only  two-dimensional  flows  will  be  considered  here, 
though  this  is  not  a strict  requirement  for  the  application  of 


1 


T 


either  method.  This  assumption  considerably  reduces  the  necessary  4 

computations  without  loss  of  any  of  the  salient  features  of  either 

approach. 

A number  of  successful  calculations  have  been  performed  which 

(3-7) 

are  capable  of  describing  transonic  flows  over  blunt  bodies 

i 

These  employ  either  the  exact  time-dependent  Euler  equations,  or 
the  full  unsteady  potential  equation  when  the  assumption  of  irro- 

tationality  is  made.  While  such  methods  are  often  time  consuming,  ^ 

they  are  able  to  properly  treat  the  time  behavior  of  bodies  oscil-  I 

J 

lating  at  high  frequency.  It  will  be  assumed  here  that  the  body  is  J 

thin  and  the  oscillation  frequency  low,  in  which  case  a limiting  1 

form  of  the  complete  potential  equation  may  be  used  to  obtain  a 1 

I 

flow  field  description.  ] 

I 

In  the  case  of  two-dimensional  inviscid  irrational  unsteady  j 

I 

flow,  the  governing  equation  for  the  velocity  potential  function  ^ 

* (f  J'"  ‘-i  : 

* |<'^  * * (0'^  (1,  ' 
3x  3x3 t 3x  3y  3x3y  3y  3y3t 

I 

If  this  equation  is  to  be  applied  in  order  to  describe  the  flow 

over  an  airfoil  with  maximum  thickness  to  chord  ratio,  6,  then  in  ' 

the  limiting  process  -»•  1,  6 -+■  0,  a much  simpler  form  results. 

A careful  study  of  the  characteristic  length  and  time  scales  for 


2 


low  frequency  oscillations  reveals  the  dominant  terms  of  Equation 
1 to  be  the  following: 

(Y  + l)(|fy/2  - [(Y 

In  actual  application,  even  more  simplification  follows  as  it  is 
only  necessary  to  include  the  dominant  nonlinear  term  in  order  to 
account  for  first  order  transonic  effects.  Equation  2 is  the  form 
considered  here  and  it  is  noted  because  the  specific  implementation 
of  one  of  the  methods  analyzed  leaves  the  time  variable  unsealed 
so  that  it  is  not  possible  to  deduce  its  exact  origin. 

For  this  study,  one  airfoil  was  selected  for  investigation  and 
was  prescribed  to  oscillate  sinusoidally  in  time  in  one  mode  of 
motion,  namely  pitch.  Flow  field  calculations  were  performed  and 
solutions  to  the  two-dimensional  unsteady  low  frequency  small  dis- 
turbance transonic  potential  equation  was  obtained  by  both  the 
method  of  harmonic  analysis  and  by  time  integration  for  several 
values  of  freestream  Mach  number  and  reduced  frequency.  Results  of 
these  calculations  are  compared.  It  should  be  noted  that  neither 
method  is  restricted  to  the  single  degree  of  freedom  (pitch)  sel- 
ected here.  The  airfoil  choosen  was  an  NACA  64A010  airfoil  which 
is  symmetric  and  ten  percent  thick,  and  is  considered  representative 
of  transonic  airfoils  currently  in  use.  Asymmetric  behavior  is  ex- 


- l)U^/2  + 


- a* 


^1! 


9$ 

9x 


9^$ 


9^<I> 

9x9t 


(2) 


= 0. 


3 


amlned  by  oscillating  the  airfoil  about  a mean  nonzero  angle  of 
attack. 
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SECTION  II 


HARMONIC  ANALYSIS 

The  method  of  harmonic  analysis  assumes  that  the  unsteady 
motion  of  an  airfoil  oscillating  with  small  amplitude  may  be  re- 
garded as  a perturbation  about  the  steady  state.  The  velocity  po- 
tential function  is  expanded  in  a series  of  Increasing  powers  of 
a small  parameter  which  is  a measure  of  the  amplitude  of  the  un- 
steady motion  of  the  boundary.  This  results  in  a sequence  of  par- 
tial differential  equations  for  the  perturbation  potentials  with 
the  zeroth  order  equation  recovering  the  steady-state  result.  All 
higher  order  equations  are  linear  with  non-constant  coefficients 
which  now  depend  upon  the  steady  solution.  This  set  of  equations 
may  be  solved  using  the  well  l<nown  relaxation  algorithms  which  have 
been  utilized  for  steady-state  transonic  flow  calculations.  Har- 
monic analysis  has  been  used  successfully  for  unsteady  airfoil  and 
rectangular  wing  motions'^  \ The  leading  order  equations  of  this 

method  are  summarized  here  for  the  specific  application  to  the 

(9) 

pitching  motion  of  an  oscillating  airfoil 
The  velocity  potential  is  written  as 

$(x,y,t)  = U„x  +u„c62/^[(l  + 4)^(5,  n.]p,  T^)  (3) 

where  c is  the  airfoil  chord  and  6 the  thickness  ratio  with 

(A) 

(5) 


5 = x/c, 

n.p  - y[(i  + y)'Sm2]*/Vc, 


5 


and 


= U^t[(l  + y)6m2]2/VcM^. 


(6) 


It  Is  assumed  that  the  angle  of  attack  a,  may  be  expressed  as 

a *>  a + a (7) 

s u ' 

where  oi  is  the  oscillation  frequency  and  a/6<<l.  The  perturbation 
potential  Is  now  decomposed  Into  steady  and  unsteady  components,  as 

n^,  T^)  = n^)  + ye^%. 

Here 


Q = kM^[(l  + 


(9) 


where 


k ■ (Jjc/U 
00 


(10) 


is  the  reduced  frequency.  It  Is  noted  that  - (Dt.  The  airfoil 
boundary  is  described  by 


y - F(x,  t)  ■ 0 for  0 < x < c 


(11) 


with 

F(x,  t)  - c5[f(0  “ aC/6],  (12) 

a given  by  Equation  7 and  f(^)  a known  function  corresponding  to 
the  airfoil  geometry. 

If  Equations  3-6  and  Equation  8 are  now  substituted  into 


6 


the  reduced  potential  Equation  2 (or  Equation  1) , the  following  ex- 
pressions for  (|)g  and  are  obtained: 


K = (1  - M2)[(l  + 

00  ^ 


(15) 


is  the  transonic  similarity  parameter.  The  corresponding  boundary 
conditions  for  Equations  13  and  14  are  obtained  from  the  flow  tan- 
gency  condition  on  the  airfoil  surface,  from  the  Kutta  condition 
at  the  trailing  edge  with  a constant  jump  in  potential  across  the 
vortex  sheet  in  the  wake,  and  by  requiring  the  perturbation  velocity 
to  vanish  far  from  the  airfoil  surface.  In  the  limit  as  6 -►  0,  the 
airfoil  and  wake  conditions  may  be  applied  on  y = 0.  Thus  the 
boundary  conditions  become 


34,s 

an. 


- a^/6  on  n. 


0 for  0 < C ^ 1» 


"1 

I 0 on  n,  - 0 for  c > 1. 


0 for  0 < 5 < 

- 0 for  5 > 1, 


(16) 

(17) 

(18) 

(19) 

(20) 
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and 


(21) 


-*•  0 as  + 


Here  the  square  brackets  denote  the  "jump"  In  the  enclosed  quantity. 

In  this  study,  for  the  purpose  of  comparison  with  time  Inte-' 
gratlon,  the  unsteady  variation  of  the  pitching  angle,  a,  was  pre- 
scribed to  have  a sinusoidal  variation  In  time;  l.e.. 


a “ a + a slnCwt) . 
s u 


(22) 


The  unsteady  flow  field  response  to  this  variation  may  be  obtained 
from  the  solutions  to  Equations  13  and  14  as 

n^,  T.J.)  - (|)g(C,  y + (a^/6)Im[4)^(5,  n^)e^^^'^].  (23) 

In  terms  of  the  perturbation  potential,  the  pressure  coefficient 
defined  In  the  usual  manner  Is  then  given  by 

IOt, 


-262/*[(1  + (a^/6)Im^^  e 


(24) 


from  which  the  corresponding  lift  coefficient  per  unit  span  may  be 
obtained : 


C.  - /I  (C  - C.)dK.  (25) 

Several  observations  concerning  this  formulation  seem  appro- 
priate here.  The  steady  small  disturbance  transonic  equation  ex- 
hibits a type  dependence  which  Is  In  Itself  a function  of  the  sol- 
ution of  the  equation;  l.e..  Equation  13  Is  hyperbolic,  parabolic, 
or  elliptic  for  ( *^8  - K ) greater  than,  equal  to,  or  less  than 


8 


zero  respectively.  Solutions  to  this  equation  have  been  efficiently 
obtained  by  relaxation  procedures  based  on  the  mixed  differencing 
scheme  of  Murman  and  Cole  (13)  which  accounts  for  the  local  nature 


of  the  flow, 
"shock  waves" 


For  a given  flow  situation,  then,  the  location  of 


will  occur  where  ( 


9(J) 

K ) changes  sign. 


In  the 


case  of  the  unsteady  perturbation  Equation  14,  type  dependence 

94) 

again  depends  only  on  the  steady  result  (i.e.,  on  - K )• 


Thus  the  unsteady  position  of  any  shock  waves  which  form  is  con- 


strained to  vary  only  slightly  from  the  steady  location.  In  fact, 
this  variation  is  proportional  to  (a^/6)  where  it  will  be  recalled 
that  a^/6<<l.  Hence,  harmonic  analysis  precludes  a description  of 
large  scale  shock  excursions  which  have  commonly  been  observed  in 
transonic  experiments.  This  results  directly  from  the  linearization 
about  the  steady  state.  Time  integration,  on  the  other  hand,  makes 
no  such  assumptions.  Comparison  of  the  results  of  these  two  meth- 
ods will  help  indicate  under  which  conditions  the  harmonic  assump- 
tions remain  valid. 

One  distinct  advantage  of  the  harmonic  method  is  that  due  to 
the  linearization,  all  aerodynamic  coefficients  now  become  linear 
functions  of  the  unsteady  deflection  angle,  ot^.  Therefore,  these 
coefficients  depend  only  on  the  reduced  frequency,  k.  Such  is  not 
true  of  the  time  integration  technique,  because  of  the  complete 


nonlinear  treatment  involved. 


SECTION  III 
TIME  INTEGRATION 

Numerical  finite  difference  schemes  that  Integrate  the  equa- 
tions of  motion  In  time  are  capable  of  simulating  nonlinear  unsteady 
transonic  flow  phenomena.  Including  Irregular  shock  wave 
motlons^^’^  6,14  17)^  These  schemes  usually  rely  upon  shock  cap- 
turing techniques  to  resolve  shock  waves,  thereby  requiring  con- 
servative differencing  form  of  the  governing  equations.  This  is 
not  always  convenient.  Time  explicit  finite  difference  procedures 
(3-6)  generally  have  a time  step  restriction  for  stability  such 
that  computations  are  often  lengthy.  More  recently,  however,  fully 
Implicit  technlques^^^  have  relieved  this  constraint,  making 
time  Integration  computationally  efficient. 

Details  of  the  form  of  the  equation  to  be  Integrated  In  time 
for  the  case  of  an  airfoil  oscillating  In  pitch  are  summarized  here. 
The  velocity  potential  Is  expanded  similarly  to  Equation  3 as 

•Kx,  y,  t)  « U^x  + U^c6*^®(()g(C,  Tig,  Tg)  (26) 

where 

- y<5‘^Vc,  (27) 

Tg  ■ a)t,  (28) 

and  C is  given  by  Equation  4.  If  these  definitions  ere  substituted 
into  Equation  2,  there  results 
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(29) 


fci  - - (1  + Y)M^  ^ '*’b+  = 

L 9?  J 8C  * 9n„2 


= 2kM 


^5-2/3 


d\ 

9?9t, 


The  airfoil  surface  is  described  by  Equations  11  and  12  such  that 
the  boundary  conditions  corresponding  to  Equation  29  are  the  fol- 
lowlne: 


9<t>B 

'^B 


df 

-r^  - a/6  on 
d^ 


0 for  0 < C < 1» 


(30) 


(31) 


(32) 


We  note  here  that  in  Equation  30,  a is  the  Instantaneous  pitch 
angle  which  is  prescribed  by  Equation  22  for  this  study.  With  this 
formulation  the  pressure  coefficient  becomes 

C - -26*  * ^'♦'b  (33) 

P 35  ’ 

and  the  lift  follows  according  to  Equation  25. 
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SECTION  IV 


DETAILS  OF  THE  NUMERICAL  SOLUTIONS 

We  now  seek  solutions  to  Equations  13  and  14  with  boundary 
conditions  16  - 21  and  Equation  29  with  boundary  conditions  30  - 
32  respectively  for  several  values  of  freestream  Mach  number,  re- 
duced frequency,  and  angle  of  attack.  Existing  computer  codes  have 
been  developed  for  the  express  purpose  of  obtaining  such  solutions 
for  the  form  of  the  equations  arrived  at  In  Sections  II  and  III. 

A certain  amount  of  freedom  Is  available  to  the  user  of  these  codes 
and  has  been  exercised  In  order  to  make  the  resultant  comparisons 
as  reliable  as  possible.  Each  code  has  certain  prominent  character- 
istics which  are  central  to  Its  respective  formulation.  These  are 
noted  In  this  section.  It  Is  the  Intention  to  achieve  a comparison 
between  the  alternative  basic  treatments  of  the  unsteady  low  fre- 
quency small  disturbance  transonic  potential  equation,  rather  than 
any  numerical  behavior  Inherent  to  a specific  computer  code.  This 
Is  not  an  easy  task  as  we  are  dealing  with  the  actual  solution  of 
finite  difference  equations  as  opposed  to  nonlinear  differential 
equations. 

Equation  13  was  solved  via  computer  code  STRANS2  (9,20)  using 
the  relaxation  procedure  mentioned  previously.  The  resulting  sol- 
ution for  (Ji  Is  then  Input  Into  code  UTRANS2  (9,20)  which  solves 
s 

Equation  14  for  (Ji^,  again  by  relaxation.  These  two  codes  were  for- 
mulated for  compatible  usage.  Time  Integration  of  Equation  29  Is 
achieved  by  employing  code  LTRAN2  (17),  which  makes  use  of  an  effl- 
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dent  time- implicit  method.  In  addition,  LTRAN2  has  the  capability 
of  solving  the  steady  transonic  equation  by  the  same  relaxation 
procedure  of  STRANS2  and  UTRANS2.  It  will  be  of  Interest  to  con- 
sider comparisons  of  some  steady  solutions. 

4.1  Airfoil  Boundary  Conditions 

STRANS2  and  UTRANS2  require  as  input  an  analytic  description 

df 

of  the  airfoil  surface  slope,  ^ . This  was  accomplished  by 

fitting  a least  squares  curve  to  NACA  tabular  airfoil  data  for  an 
assumed  functional  dependence  of  f on  ? (21) . The  function  was 
selected  as 

f(0  = b ^ + i b (34) 

n=0  ^ 

1/2 

from  which  the  slope  follows  directly.  We  note  the  term  C is 
Included  to  approximate  the  blunt  leading  edge,  thus,  making  the 
surface  slope  infinite  at  C = 0.  Small  disturbance  theory  is  in- 
capable of  treating  this  singularity  and  in  general  either  the 
leading  edge  is  not  Included  in  the  computational  mesh,  or  the 
leading  edge  slope  is  modified  to  some  finite  value  in  the  limit 
C ->  0,  in  order  to  accomodate  blunt  nosed  airfoils.  The  former 
technique  was  employed  here  for  both  harmonic  analysis  and  time  in- 
tegration. LTRAN2  approximates  the  airfoil  slope  by  passing  a cubic 
spline  through  the  tabular  data  and  extracting  the  corresponding 
slope. 

4.2  Far field  Boundary  Conditions 

The  harmonic  farfleld  boundary  condition.  Equations  18  and  21, 
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are  treated  by  obtaining  the  limiting  form  of  the  potentials  (iji^ 
and  (J>^)  far  from  the  airfoil  and  applying  this  result  at  a large 
but  finite  distance.  This  approach  was  used  in  STRANS2  and  UTRANS2 
with  considerable  savings  in  the  extt  " of  the  computational  mesh. 
For  the  steady  potential,  ij)^,  the  limiting  form  is  the  farfield 
condition  derived  by  Klunker  (22),  which  results  from  application 
of  Green's  theorem  to  Equation  13.  In  a similar  fashion,  a far- 
field expression  for  (J)^  was  obtained  in  terms  of  the  Hankel  func- 
tion of  the  second  kind  of  order  zero  (9).  A finite  series  expan- 
sion is  used  to  approximate  this  Hankel  function  in  UTRANS2.  If  the 
farfield  condition  is  applied  too  far  from  the  body  it  was  found 
that  this  series  may  fall  to  converge.  This  fact  was  of  Importance 
when  numerical  meshes  were  chosen  as  will  be  shown  in  the  next  sub- 
section. 

For  time  integration,  no  approximate  farfield  condition  has 
been  found.  Therefore,  LTRAN2  applies  the  mean  steady-state  lim- 
iting form  of  the  potential  extremely  far  from  the  body.  While 
this  treatment  is  not  exact,  it  has  produced  no  numerical  anomalies. 
The  farfield  boundary,  however,  must  be  taken  sufficiently  far  from 
the  airfoil  such  that  no  waves  reflected  from  this  boundary  reach 
the  immediate  vicinity  of  the  airfoil  over  the  interval  of  time 
covered  by  the  integration. 

4.3  Computational  Meshes 

Codes  STRANS2  and  UTRANS2  have  the  capability  of  accepting  in- 
put data  for  both  Che  number  and  distribution  of  numerical  mesh 
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points.  This  is  not  in  general  true  of  LTRAN2.  As  the  spatial 
differencing  of  the  finite  difference  equations  of  the  time  inte- 
gration in  LTRAN2  is  similar  to  that  employed  by  the  harmonic  anal- 
ysis codes,  then  if  the  computational  meshes  are  identical,  the 
truncation  error  incurred  in  the  respective  solutions  should  be 
formally  of  the  same  order.  Because  the  scaling  of  the  y-coordinate 
is  not  the  same  for  both  methods  (see  Equations  5 and  27) , the 
meshes  were  made  identical  by  forcing  values  of  n„  and  n,,  to  cor- 
respond  to  the  same  value  of  the  physical  variable,  y.  Identical 
values  of  C were  used  in  all  programs. 

The  meshes  employed  in  STRANS2  and  UTRANS2  were  truncated  ver- 
sions of  the  one  used  in  LTRAN2  because  of  the  treatment  of  the  far- 
fleld  boundary  conditions  described  in  the  last  subsection.  Both 
meshes  used  33  points  distributed  along  the  airfoil  surface  with  a 
point  at  the  trailing  edge  and  the  leading  edge  situated  between 
mesh  points.  The  minimum  C”Step,  A C min,  was  0.00330  at  the  lead- 
ing edge.  Minimum  mesh  sizes  in  the  normal  direction  occured  at 
y “ 0 with  a symmetrical  distribution  of  points  above  and  below 
this  axis.  Mesh  boundaries  and  other  pertinent  data  are  given  be- 
low: 

a.  STRANS2,  UTRANS2 

(1)  number  of  C-points  ■ 70 

(2)  number  of  n.p-polnt8  ■ 43 

^3)  - 0-02151 

(4)  -8.19140  < P < 9.02411 
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(5)  -10.13300  1 1 10.13300 

b.  LTRAN2 

(1)  number  of  C-polnts  = 99 

(2)  number  of  n„-points  = 79 

D 

(3)  = 0.02000 

(4)  -1033.53047  < C 1 855.91313 

(5)  -811.12200  < Ti„  < 811.12200 
4.4  Accuracy  of  the  Calculations 

Steady-state  solutions  for  provided  by  STRANS2  were  con- 
sidered converged  when  the  variation  In  41^  at  all  mesh  points  be- 
tween consecutive  Iterations  was  less  than  10  This  same  criteria 
was  applied  when  steady  solutions  were  obtained  from  LTRAN2.  In 
the  case  of  the  complex  potential,  (Ji^,  a less  severe  standard  was 
Invoked.  This  was  because  of  the  extremely  slow  rate  of  convergence 
of  the  solution  procedure  employed  by  UTRANS2  for  solving  Equation 
14.  More  accuracy  In  4)^  was  Insisted  upon  because  of  Its  appear- 
ance In  Equation  14.  Unsteady  solutions  were  accepted  as  converged 
when  the  variation  In  1 41  I between  consecutive  Iterations  was  less 

I Ty  I 

-4  I I 

than  10  . Here,  |4i^|  Implies  the  absolute  value  In  a complex 

sense.  This  was  considered  more  than  sufficient  for  an  adequate 
comparison  with  time  Integrated  results. 

The  steady-state  solutions  obtained  from  LTRAN2  were  used  as 
Initial  profiles  for  all  time  Integrations.  Starting  with  these 
values,  solutions  were  advanced  In  time  Invoking  the  forced  boundary 
condition.  Equation  22.  The  time  step  size  used  for  any  given  case 


16 


corresponded  to  360  time  steps  over  one  period  of  forced  oscillation 
at  the  reduced  frequency.  It  was  determined  that  such  a time  step 
was  sufficiently  small  by  increasing  this  value  by  up  to  a factor 
of  five  without  variation  in  the  first  three  significant  figures 
of  the  computed  lift  and  local  pressure  coefficients. 

After  solutions  were  advanced  sufficiently  in  time,  the  effect 
of  the  initial  conditions  became  negligible  and  the  solutions 
attained  a periodic  behavior  due  to  the  periodic  boundary  condition. 
It  is  this  periodic  behavior  which  we  desire  to  compare,  not  the 
small  time  history  which  can  never  be  predicted  by  harmonic  analysis. 
Solutions  were  considered  periodic  when  repeatability  to  three  sig- 
nificant figures  was  attained  over  two  successive  periods  by  all 
local  pressure  coefficients.  For  the  cases  examined  in  this  study, 
effects  of  the  initial  conditions  usually  persisted  for  less  than 
one  period.  As  a result  all  solutions  were  forced  for  three  periods 
of  oscillation  starting  from  the  initial  profile. 

Lift  coefficients  for  both  methods  were  computed  according  to 
Equation  25  by  numerical  integration.  Codes  STRANS2  and  UTRANS2 
employed  trapezoidal  rule  for  this  purpose,  whereas  LTRAN2  used  a 
Simpson's  rule.  Although  this  renders  the  LTPAN2  results  more  accu- 
rate, both  techniques  were  considered  sufficiently  exact  for  the 
comparisons  made  here. 
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SECTION  V ^ 

RESULTS  ' 

1 

Figure  1 shows  the  comparison  of  the  steady-state  results 
from  STRANS2  and  LTRAN2  in  terms  of  the  surface  pressure  coeffi- 
cient for  the  case  M =0.72  and  a = 0.  Because  the  airfoil  is 

00 

symmc'tric , the  pressure  distributions  on  the  upper  and  lower  sur- 
faces are  identical.  It  is  noted  that  this  case  is  subcritical. 

\greement  between  the  two  results  is  seen  to  be  quite  good  over  the 

last  50%  of  the  airfoil.  Lack  of  agreement  on  the  front  portion  ' 

J 

of  the  surface  is  attributed  to  the  alternative  methods  used  to 
prescribe  the  airfoil  slope  by  the  respective  codes.  In  particular, 
it  is  believed  that  the  anomalous  behavior  near  C “ 0-1  in  tha 
LTRAN2  solution  results  from  the  cubic  spline  used  to  approximate 
the  airfoil  surface.  Other  than  this  behavior,  the  results  compare 
reasonably  well. 

The  unsteady  lift  coefficient  as  a function  of  time  for  a 1“ 
pitching  oscillation  at  a reduced  frequency  of  0.05  and  freestream 
Mach  number  of  0.72  is  displayed  in  Figure  2.  Here  the  time  origin 
is  purely  fictitious  and  has  been  Introduced  solely  for  the  purpose 
of  comparison.  Respective  time  histories  of  lift  calculated  by  each 
method  were  aligned  at  the  maximum  and  minimum  lift  values  for  this, 
and  subsequent  figures.  The  results  of  the  two  methods  are  seen  to 
compare  exceptionally  well  except  near  the  peak  lift  value.  Even 
there,  the  predicted  results  are  with  5%  of  each  other.  A compar- 
ison of  the  unsteady  surface  pressure  coefficients  corresponding 
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o/t  (DEGREES) 


Figure  2.  Comparison  of  Unsteady  Lift  Coefficient 

for  a 1*  Pitching  Oscillation  at  = 0.72 
and  k = 0.05. 
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to  this  case  are  shown  in  Figure  3.  The  time  here  was  chosen  as 
, that  at  which  the  maximum  lift  occured  where  the  widest  variation 

between  the  two  results  appeared.  It  Is  noted  that  the  flow  has 
remained  subcrltlcal.  Generally  the  results  compare  quite  favorably 
with  the  greatest  discrepancy  between  the  two  methods  agaln^  occur- 
ring over  the  front  portion  of  the  surface.  The  anomalous  behavior 
In  the  Initial  profile  (steady-state)  of  the  time  Integration  appears 
to  have  persisted  Into  the  unsteady  result. 

In  Figure  4 the  unsteady  lift  coefficient  at  the  same  Mach 
number  and  pitching  displacement  Is  displayed  for  a higher  reduced 
frequency,  k ° 0.20.  This  comparison  Is  quite  similar  to  that  seen 
at  the  lower  frequency  (Figure  2)  with  the  maximum  disagreement  be- 
tween the  two  results  occurring  at  the  peak  lift  value.  It  Is  In- 
teresting to  note  that  both  harmonic  analysis  and  time  Integration 
Indicate  lower  values  for  the  lift  coefficient  at  this  higher  fre- 
quency (k  « 0.20)  than  were  predicted  for  k * 0.05.  In  addition, 
for  k • 0.20  time  Integration  predicts  a lower  lift  than  harmonic 
analysis.  The  opposite  behavior  was  observed  for  k - 0.05  (Figure 
2).  Corresponding  unsteady  surface  pressure  coefficients  for  this 
case  are  shown  In  Figure  5.  Again  these  results  are  quite  similar 
to  those  of  the  lower  frequency  case  shown  In  Figure  3. 

As  a representative  example  of  a comparison  of  supercritical 
calculations  by  the  two  computational  methods,  a freestream  Mach 
number  of  0.82  was  selected.  Steady-state  surface  pressure  coeffi- 
cients are  shown  In  Figure  6.  A well-defined  shock  Is  seen  to  lie 
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between  C ~ 0.5  and  ^ <■  0.6.  The  same  anomalous  behavior  of  the 
LTRAN2  result  for  the  case  ■ 0.72  is  present  here.  Downstream 
of  this,  a slight  variation  in  the  profiles  up  to  the  mid-chord 
is  noted,  with  good  agreement  from  there  rearward. 

Figure  7 Indicates  the  unsteady  comparison  at  this  Mach  number 
for  a 1®  pitching  oscillation  at  a reduced  frequency  of  0.05  in 
terms  of  the  lift  coefficient.  Unlike  the  lower  Mach  number  results, 
there  is  now  considerable  difference  in  the  two  solutions.  At  the 
peak  lift  value  the  time  integration  is  seen  to  produce  approximately 
50%  more  lift  than  the  harmonic  analysis.  Clearly,  this  effect  is 
due  to  the  presence  of  the  shock  and  its  associated  motion. 

The  unsteady  surface  pressure  coefficients  at  the  value  of  time 
corresponding  to  the  maximum  lift  for  this  supercritical  case  are 
displayed  in  Figure  8.  Near  C “ 0.1,  the  anomalous  behavior  in  the 
time  integration  appears  to  have  become  more  pronounced.  Time  in- 
tegration has  predicted  a lower  and  a higher  maximum  pressure  on 
the  lower  surface  than  harmonic  analysis.  In  addition,  because  the 
airfoil  is  symmetric,  the  average  of  the  position  of  the  shock  on 
the  upper  and  lower  surfaces  corresponds  to  the  time  mean  shock 
location.  This  is  somewhat  further  aft  for  the  time  integration 
than  for  harmonic  analysis  which  always  predicts  the  mean  location 
to  coincide  with  the  steady-state  value.  Others  have  observed  this 
identical  behavior  (23).  It  is  also  evident  that  harmonic  analysis 
has  constrained  the  shock  motion  from  moving  as  far  aft  on  the 
upper  surface  or  as  far  forward  on  the  lower  surface  than  is  permitted. 
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Figure  8.  Comparison  of  Unsteady  Surface  Pressure  Coefficient  at  Time  Corresponding 
to  Maximum  Lift  for  a 1“  Pitching  Oscillation  at  M = 0.82  and  k = 0.05. 


by  integrating  in  time. 

In  order  to  examine  asymmetric  behavior  of  the  solutions  gen- 
erated by  the  two  methods,  the  airfoil  was  oscillated  about  a non- 
zero angle  of  attack.  Steady-state  surface  pressure  coefficients 
at  = 0.72  and  a = 1“  are  displayed  in  Figure  9.  This  case  is 
seen  to  be  sub-critical  and  it  is  interesting  to  note  that  this 
steady  result  is  quite  similar  to  the  two  unsteady  results  at  the 
same  Mach  number  (Figures  3 and  5).  The  time  varying  lift  coeffi- 
cient for  a 0.5“  pitching  oscillation  about  the  steady  1“  angle  of 
attack  for  a reduced  frequency,  k = 0.05,  is  shown  in  Figure  10. 
Here  the  shape  of  the  time  history  lift  curve  generated  by  each 
method  is  virtually  identical.  However,  time  integration  predicts 
a higher  time  mean  lift  coefficient  than  that  Indicated  by  harmonic 
analysis.  This  corresponds  to  the  same  behavior  as  was  observed 
for  the  case  = 0.82  where  it  was  noted  that  harmonic  analysis 
always  recovers  steady-state  conditions  as  the  time  mean  values, 
which  is  not  necessarily  true  of  time  integration.  Figure  11  in- 
dicates the  comparison  of  surface  pressures  for  this  case  at  the 
time  corresponding  to  maximum  lift.  The  flow  has  remained  just 
subcrltlcal,  with  most  of  the  discrepancy  between  the  two  results 
occurring  near  the  leading  edge.  Abrupt  changes  in  the  slope  of 
the  pressure  curve  which  are  evident  in  time  Integrated  profile  may 
have  resulted  in  loss  of  accuracy  in  the  calculation  of  the  lift 
coefficient  through  numerical  integration. 
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Figure  10.  Comparison  of  Unsteady  Lift  Coefficient  for  a 
0.5®  Pitching  Oscillation  about  a Mean  Angle 
of  Attack  of  1"  at  M_  = 0.72  and  k = 0.05. 


SECTION  VI 


CONCLUSIONS 

Solutions  to  the  unsteady  low  frequency  small  disturbance 
transonic  potential  equation  have  been  obtained  by  both  the  method 
of  harmonic  analysis  and  by  the  method  of  time  Integration  for  the 
pitching  oscillation  of  an  NACA  64A010  airfoil  at  several  combin- 
ations of  freestream  Mach  number,  reduced  frequency,  and  displace- 
ment amplitude.  As  far  as  was  possible,  computational  differences 
In  obtaining  those  solutions  were  minimized.  When  this  could  not 
be  accomplished,  these  differences  were  noted.  It  Is  considered 
significant  that  all  solutions  were  generated  using  an  Identical 
computational  mesh  In  the  Immediate  vicinity  of  the  airfoil.  The 
choice  of  free  parameters  corresponding  to  each  particular  case 
considered  (M^,  k,  pitching  amplitude)  was  specifically  restricted 
so  as  not  to  violate  the  Inherent  assumption  of  either  method.  The 
reduced  frequency  was  never  more  than  0.20  as  the  low  frequency 
approximation  becomes  a poor  assumption  above  a value  of  about  0.40. 
Pitching  amplitudes  were  specified  to  be  1®  or  less  so  as  not  to 
violate  the  small  disturbance  assumption.  While  harmonic  analysis 
was  not  expected  to  predict  large  scale  shock-wave  motions  with 
accuracy,  a supercritical  case  was  considered  In  order  to  deduce 
the  range  of  acceptability  for  this  method  of  analysis. 

Comparisons  between  the  solution  methods  were  provided  In  terms 
of  the  time  dependent  lift  coefficient  and  the  unsteady  surface 
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pressure  coefficient.  For  oscillations  in  a subcritical  flow  field 
, about  a zero  mean  angle  of  attack  the  two  methods  were  found  to 

j agree  quite  well.  When  a large  shock  appeared,  harmonic  analysis 

failed  to  reproduce  the  unsteady  shock  motion  predicted  by  time  in- 
. tegration.  In  the  case  of  oscillation  about  a non-zero  mean  angle 

of  attack,  the  unsteady  lift  obtained  from  the  harmonic  result  varied 

by  a simple  translation  from  that  resulting  from  Integrating  in 
time. 
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